Convergent resummed linear 5 expansion in the critical O(N) {(j)f)\d model 
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The nonperturbative linear S expansion (LDE) method is applied to the critical O(N) cj> 4 three- 
dimensional field theory which has been widely used to study the critical temperature of condensation 
of dilute weakly interacting homogeneous Bose gases. We study the higher order convergence of the 
LDE as it is usually applied to this problem. We show how to improve both, the large- N and finite 
N — 2, LDE results with an efficient resummation technique which accelerates convergence. In the 
large N limit, it reproduces the known exact result within numerical integration accuracy. In the 
finite N = 2 case, our improved results support the recent numerical Monte Carlo estimates for the 
critical transition temperature of Bose-Einstein condensation. 
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Perturbative methods are well established tools to 
study diverse aspects of physical systems. However, in 
many important problems one is not allowed to use per- 
turbation theory, or it just breaks down. In field theory, a 
well known example concerns the description of physical 
phenomena close to critical points of continuous or sec- 
ond order phase transitions. In this situation the diver- 
gence of length and time scales of the fluctuations, associ- 
ated to infrared (IR) divergences and critical slow-down, 
respectively, result in the singular behavior of many phys- 
ical quantities, like correlation lengths, susceptibilities 
and many other parameters. In these cases one must 
recur to nonperturbative methods, such as the renormal- 
ization group methods, the e-expansion, the 1/N approx- 
imation and other techniques (for a review, see Ref. 0). 

A timely important problem associated to the break- 
down of perturbation theory near a critical point is the 
evaluation of the critical transition temperature for a 
weakly interacting homogeneous Bose gas. This appar- 
ently simple, but highly nontrivial task has been the 
source of controversy for many years. Recently, its func- 
tional form was found to be T c = Tq{1 + cian 1 / 3 + 
[c' 2 lr^an 1 / 3 ) + c 2 ']a 2 n 2 / 3 + 0(a 3 n)}, where T is the ideal 
gas condensation temperature, a is the s-wave scatter- 
ing length, n is the density and Ci,C2,c 2 ' are numerical 
coefficients However, a strong debate concerns the 
values of the numerical coefficients, especially ci, which 
has been computed by different authors, using different 
methods. Some analytical predictions included the self- 
consistent resummation schemes (ci ~ 2.90) ||, the 1/iV 
expansion at leading order (ci ~ 2.33) (i) and at next 
to leading order (ci ~ 1.71) || and also the linear 5- 
expansion (LDE) || fj| at second order (ci ~ 3.06) ||. 
The numerical methods include essentially Monte Carlo 
lattice simulations (MC). The most recent MC results are 
reported by the authors of Ref. @ (c x = 1.29 ±0.05) and 
of Ref. f0j (ci = 1.32 ±0.02). The problem is that these 
coefficients (except c 2 ) are sensitive to the IR physics at 
the critical point and so, no perturbative approach can 
be used to predict them. 

At first, one may believe that, due to the complex na- 



ture of the problem, it is unlikely that a definitive an- 
alytical prediction of the IR coefficients could be made. 
In general, higher order computations quickly become 
prohibitively difficult with the traditional nonperturba- 
tive analytical methods which, like the 1/N expansion, 
rely on the resummation of an infinite number of partic- 
ular contributions. At the same time, the LDE great ad- 
vantage is the fact that all calculations, including renor- 
malization, are performed in a perturbative way. This 
means that one always deals with a finite set of con- 
tributions at each order. This advantage is well illus- 
trated in a recent work by some of the present authors 
]ll[ , where the LDE was extended to order-<5 4 . The 
results give strong indications that, when properly ap- 
plied, the LDE leads to very precise analytical predictions 
for the nonperturbative coefficients c\ and c' 2 '. In fact, 
the results of Ref. [O for these coefficients, at second 
(ci ~ 3.06, 4 ~ 101.4), third (ci ~ 2.45, d{ ~ 98.2) and 
fourth orders (ci ~ 1.48, c' 2 ' ~ 82.9) seem to, roughly, 
converge to the precise MC results of Ref. M, where 
c' 2 ' = 75.7 ± 0.4 was also predicted. 

The important question regarding the LDE conver- 
gence properties has been addressed in the context of the 
anharmonic oscillator (AO) at zero temperature where 
rigorous proofs have been produced BL Those proofs 
have been extended to the finite temperature domain by 
Duncan and Jones |l3| who considered the AO partition 
function. Very recently, the convergence proof has been 
extended to renormalizable quantum field theories at zero 
temperature ]l4| |. However, it would be interesting to 
probe convergence in the vicinity of a phase transition 
such as for the Bose-Einstein condensation (BEC) model 
considered here. In most applications one can establish 
simple relations in between the LDE and other nonper- 
turbative methods already at order-J where one-loop di- 
agrams are present. In fact, one can show that in those 
cases the LDE either exactly reproduces 1/N results or 
produce very close numerical estimates |l2|, [ll| . Here, the 
BEC problem poses an additional difficulty since the first 
non-trivial contributions start at the two-loop level in the 
self-energies. This is a consequence of the Hugenholtz- 
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Pines theorem which eliminates the one-loop momentum 
independent contributions. In this case, it is not easy to 
establish simple analytical relations like those given in 
Refs. [|l2[ |l5| and the problem must be treated differ- 
ently. 

Braaten and Radescu in Ref. jl6| have recently re- 
visited the LDE application to the BEC problem also 
considering the convergence problem. One of the differ- 
ences between our present work and also Refs. [||. Q 
from their approach regards the choice of physical quan- 
tity to be extremized. Moreover, the coefficient c'2 and 
convergence acceleration, considered by us, were not in- 
vestigated in Ref. [jl6j. We discuss both approaches, in 
more details, in a longer companion paper |l7]] . 

In this Letter, our aim is to mainly investigate the con- 
vergence of the LDE, in the critical O(N) {4> 2 )\ d model, 
by considering its behavior at the large-iV as well as at 
the finite N = 2 limits. This study will be combined 
with a simple, but powerful, all order LDE resummation 
technique [ 18j so that we can explicitly show convergence. 

At the critical point, one can describe a weakly in- 
teracting dilute homogeneous Bose gas by an effective 
action analogous to a 0(2) scalar field model in three- 
dimensions given by 
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where <p is a two-component real scalar field. The pa- 
rameters r and u are related to the original parameters 
of the nonrelativistic action, the chemical potential /i, 
a, atomic mass m and temperature T by r — —2m/j, 
and u = A&namT ||, ||. The leading order coefficient 
of the critical temperature shift can be expressed as 
§ a = -1287T 3 [C(3/2)]- 4 / 3 A(0 2 )/u, where A<0 2 } = 
(<p )« — {4> 2 )q- The subscripts u and mean that the 
field fluctuations are to be evaluated in the presence and 
in the absence of interactions respectively. 

The implementation of LDE within this model is re- 
viewed in previous applications H (see also Ref. jl9[). 
In practice, one considers the original theory Eq. (If) 
adding a quadratic (in the fields) term 1/2(1 — <5)?7 2 </r 
where ry is an arbitrary mass parameter. At the same 
time u — ► 6u and r — * dr. One ends up with an interpo- 
lated theory described by the propagator (q 2 + rj 2 )^ 1 and 
vertices —5u , Srj 2 and — Sr. A physical quantity <E>( fe ) is 
then evaluated to an order-<5 fc using these new Feynman 
rules and following the program, including renormaliza- 
tion, of ordinary perturbation theory. This method op- 
erates for any TV and is free from IR divergences. Non- 
perturbative results for the original (^-independent) the- 
ory are obtained by optimizing with respect to rj at 
5=1. One possible optimization procedure, adopted in 
a part of this work, is the Principle of Minimal Sensitivity 
(PMS) which requires d^ k) /dr] = @. For comparison 
and to check the consistency of the optimization proce- 
dure for reproducing physically sensible results, we also 
consider the alternative optimization scheme known as 



the Fastest Apparent Convergence criterion (FAC) where 
one requires that, at a given order k, the fc-th coefficient 
of the perturbative expansion be zero. As we will see, 
both methods lead to equivalent results. 

Let us first consider, in the large-A^ limit, the usual 
analytical way ||, [ll] to define a non-trivial LDE, by 
performing an order by order perturbative evaluation 
of ((/) 2 )^ which, moreover, has the advantage of being 
straightforwardly generalized for arbitrary N values. The 
convergence property issues can be studied by comparing 
the results with the "exact" 1/N result f| c x = 2.328. 
The details of such a calculation follow from the meth- 
ods used in Ref. |p"i"|| and are given in the accompanying 
paper [|l7|. One obtains 
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where the coefficients are given by C\ = 7.249 x 10 5 , 
C* 2 = 2.032 x 10~ 6 , C 3 = 6.400 x 10" 8 , C 4 = 2.080 x 10~ 9 , 
C* 5 = 5.021 x 10~ n , C 6 = 2.760 x 10" 12 , C 7 = 3.580 x 
1CT 14 , C 8 = 6.500 x 10- 16 , C 9 = 1.090 x 10~ 17 , C 10 = 
1.040 x 10~ 19 , C n = 7.030 x 10~ 22 , C 12 = 2.810 x 10~ 24 , 
Ci 3 = 7.300 x 10~ 26 , Cu = 2.800 x lO" 28 , C 15 = 
6.000 x IO- 31 , Ci 6 = 5.780 x 10~ 33 , C 17 = 1.100 x lO 35 , 
Ci 8 = 1.130 x 10~ 37 and C 19 = 1.390 x 10~ 40 . All 
coefficients for i > 2 in Eq. (^) were obtained from i- 
dimensional integrals over Feynman parameters, that we 
have performed by using the well-known Monte Carlo 
multidimensional integration routine VEGAS pl]| . We 
should note that for such high dimension and compli- 
cated integrals, the Monte-Carlo statistical integration 
error cannot be reduced well below the percent level and 
this will be reflected on our final results. 

Here, as in the original LDE applications to the BEC 
problem §, 0, our strategy is to evaluate ((j> 2 )^ per- 
turbatively to order-<5 fe and then to extremize this quan- 
tity following the PMS. By setting u = in the optimal 



)u one immediately obtains the optimal 



and, 



at the same time, the optimal A((j) 2 )^ k K By applying the 
PMS to this quantity, at each order, one obtains solu- 
tions that can be grouped into complex families whose 
first member is real as discussed in the anharmonic os- 
cillator analogous studies (see Belief, Garcia and Neveu 
in Refs. |||) and in Ref. 0. To order-<5 20 , our results 
are shown in Table I, together with those obtained by 
the FAC optimization procedure, for the best converg- 
ing family of solutions (see Ref. jlTj for details) . In the 
very last line of Table I, we indicate the corresponding 
statistical integration accuracy as provided by VEGAS. 
Note that both, the PMS and FAC, results converge to 
similar values, attesting that our results are not merely 
an artifact of the optimization procedure. 

As far as the predictions of Ref. Jl(| are concerned we 
note that our results, optimized via PMS or FAC, are 
also very stable. Like those authors, at about 18 th order 
we achieve ~ 7% accuracy. Let us now present a way 
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TABLE I: PMS, FAC and CIRT-PMS results for ci at large- 
N, at different orders k, obtained with the best converging 
families (real part) toward the exact result ci = 2.328. 



k 


PMS 


FAC 


CIRT-PMS 


2 


4.326 


4.996 


_ 


3 


3.760 


3.774 


_ 


4 


4.400 


3.298 


4.886 


5 


3.260 


3.126 


4.258 


6 


3.000 


2.940 


3.622 


7 


2.852 


2.810 


2.450 


8 


2.760 


2.738 


2.310 


9 


2.700 


2.684 


2.340 


10 


2.652 


2.638 


2.346 


11 


2.600 


2.604 


2.346 


12 


2.580 


2.574 


2.346 


13 


2.540 


2.552 


2.346 


14 


2.538 


2.532 


2.346 


15 


2.520 


2.516 


2.346 


16 


2.504 


2.504 


2.346 


17 


2.488 


2.492 


2.346 


18 


2.480 


2.482 


2.346 


19 


2.476 


2.474 


2.346 


20 


2.468 ±0.025 


2.466 ± 0.025 


2.346 ±0.023 



of improving the above results with an efficient LDE re- 
summation technique. Performing the usual LDE inter- 
polation with -q* = 77(1 — (5) 1 / 2 and u — > Su, expanded to 
order p, defines a partial sum ^ p ^(rj, u, S) = X)n=o s ™^"> 
which for 8 — ► 1 is given formally, from the simple pole 
residues, as 



Idvexp(v) v a = 2ni/T(~a) , (5) 

one sees that the main effect of this resummation is to 
divide the original expansion coefficients at order k by 
terms F(l±fc/2) ~ (fc/2)! for large k. This damping 
of the perturbative coefficients at large order, as implied 
by this specific resummation, is rather generic and was 
exploited recently in the completely different context of 
asymptotically free models |L4| where it was shown to 
accelerate convergence of the LDE. When applied to the 
anharmonic oscillator, it is in fact (asymptotically) equiv- 
alent to the more direct LDE resummation with an order- 
dependent rescaling of the arbitrary mass, as employed 
in some of the Refs. jl2| to establish rigorous convergence 
of the LDE for the oscillator energy levels, which is it- 
self an extension of the order-dependent mapping (ODM) 
resummation technique . The contour integral resum- 
mation is however very convenient since, algebraically, it 
is simpler than the direct LDE summation, in particular 
to recover the original theory for 77* — ^ 0. The results ob- 
tained through this contour integral resummation tech- 
nique (CIRT) applied to Eq. (§ are also shown in Table 
I. Fast convergence happens already at order-^ 10 , within 
~ 1% of the exact large- N result for c\. Note the com- 
plete stability of results from this order onwards. Within 
the intrinsic numerical integration error involved in the 
computation of the coefficients in Eq. (^) we can con- 
clude on the actual convergence of the CIRT towards the 
exact 1/N result. Note how the CIRT-PMS predictions 
quickly become more accurate, than the ordinary PMS 
(FAC) results, when more approximants are considered. 

We finally turn to the finite N = 2 case, for which the 
quantity ((j) 2 )^ has been evaluated, up to order-(S 4 , in 
Ref. [jll) . Its perturbative expansion is 
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where the anticlockwise contour encircles the origin. 
Now, one performs a change of variables for the rel- 
evant 5 — ► 1 limit: 5 = 1 — v /p, together with a similarly 
order-dependent rescaling of the arbitrary mass parame- 
ter, 77 > 77 p 1 / 2 , where the power 1/2 is dictated by the 
scalar mass interaction term. For p — > 00 this resumma- 



tion can be summarized as the replacement 77* 
followed by the contour integration 
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where the coefficients are given by K\ — 3.222 x 10 5 , 
K 2 = 1.524 x 1CT 6 and K 3 = 1.042 x 10~ 7 . The 
contour integral technique applied at large-iV also im- 
proves the finite N = 2 result of Ref. jll|] showing 
that the complex family of solutions has real parts that 
converge to ci = 1.15 ± 0.01, whereas the most recent 
Monte Carlo predictions are || c% = 1.29 ± 0.05 and 
p| a = 1.32 ± 0.02. Apart from {<p 2 ) {k \ the quantity 
<5rl fe ) _ — Sren(0) also enters the evaluation of the order- 
a 2 coefficient c 2 ' which appears in the T c expansion [0. 
For N = 2 its order-5 4 perturbative evaluation is 0] 



where the "weight" exp(v)/v originates from dS (I — 
S)^ 1 — > dv/v; lim p ^ 00 (l — v/p)^ 13 ^ 1 — exp(Tj), and the 
original contour was deformed to encircle the branch cut 
Re[v] < 0. Noting that the LDE produces a power series 
in ti( fe+1 ) /r] k and by using 



SrW = 5 ^L + S 2 u 2 A 2 
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M 



0.59775 



,3 u3 a. , xi u4 



6*—A 3 + 5 4 



A 4 + 0(S 5 ) , 



(7) 
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where the coefficients are A2 — 1.407 x 10 3 , A3 = 
8.509 x 10~ 5 and A 4 = 3.523 x 10~ 6 . Treating Eq. 

(0) with the CIRT one obtains the result Re [r^ 4) (M = 
u/3)} = 0.0010034u 2 which, together with the CIRT im- 
proved (0 2 }1 4) result, leads to d{ = 84.9 ± 0.8 |7), while 
the Monte Carlo estimate is 0% = 75.7 ± 0.4 [2j . Note 
that the scale M = u/3 was originally chosen in those 
Monte Carlo applications. 

In conclusion, the LDE has been applied successfully to 
many different problems in field theory where standard 
perturbation theory does not apply. But despite of its 
successes, its applicability to higher orders and the study 
of its convergence properties in field theory, beyond the 
anharmonic oscillator problem ]l2[| , have proven to be a 
challenge. Here, we have studied the convergence of the 
LDE as applied to the critical point of Bose-Einstein con- 
densation. We have shown that the perturbative LDE for 
the large- iV case converges towards the exact result, once 
resummed to all orders. We have also used the LDE to 
explicitly evaluate, in the large-iV limit, the coefficient 
ci to order-5 20 using two different optimization proce- 
dures (PMS and FAC), leading to equivalent converging 
results, but in both cases ~ 7% accuracy is only achieved 
at about order-(5 18 . Then, we have shown how the pow- 
erful contour integral resummation technique (CIRT) ac- 
celerates convergence already at order-<5 10 within 1% ac- 
curacy or less, which appears in fact only limited by the 
intrinsic accuracy of the Monte-Carlo integration used to 
evaluate the coefficients of the relevant series. This same 



technique was extended to the relevant finite N = 2 case 
where the recent order-<5 4 results for c\ and cJ,' |llj] were 
improved. It is worth remarking that, the c\ values 3.06, 
2.45 and 1.48 obtained at orders-(5 2 , S 3 and S 4 in Ref. 
Jll| are close to 2.90, 2.33 and 1.71 obtained by rcsum- 
ming "setting sun" contributions Q , leading order || and 
next to leading order 1/N || contributions, respectively. 
The similarity between the LDE values at a given order 
and the values produced by each one of those analytical 
nonperturbative approximations should come as no sur- 
prise if one considers the type of graphs resummed by 
the LDE, at each order. At order-5 2 , only "setting sun" 
contributions are considered while more typical 1 /N lead- 
ing order contributions appear at order-5 3 . At order-<5 4 , 
graphs which would belong to the 1/N at next to leading 
order also contribute. This simple consideration shows 
that the hierarchy of LDE numerical values is not a mere 
coincidence but a consequence of the type of graphs con- 
sidered at each order. The numerical differences are due 
to the fact that the LDE actually mixes up those contri- 
butions, at each order in <5, irrespective of their 1/N or- 
der. Here, apart from supporting the recent Monte Carlo 
predictions, our improved results reinforce the potential 
of the LDE as a powerful tool to treat nonperturbative 
problems in field theory. 
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